Thermodynamics of magnetized binary compact objects 



O 

(N 

^— > 
O 

O 

(N 



Koji Uryu, 1 Eric Gourgoulhon, 2 and Charalampos Markakis 3 

Department of Physics, University of the Ryukyus, Senbaru, Nishihara, Okinawa 903-0213, Japan 
2 Laboratoire Univers et Theories, UMR 8102 du CNRS, 
Observatoire de Paris, Universite Pans Diderot, F-92190 Meudon, France 
3 Department of Physics, University of Wisconsin- Milwaukee, P.O. Box 413, Milwaukee, WI 53201 

(Dated: 28 August 2010) 

Binary systems of compact objects with electromagnetic field are modeled by helically symmet- 
ric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived 
thermodynamic laws for helically-symmetric perfect-fluid spacetimes are extended to include the 
electromagnetic fields, and electric currents and charges; the first law is written as a relation be- 
tween the change in the asymptotic Noether charge SQ and the changes in the area and electric 
charge of black holes, and in the vorticity, baryon rest mass, entropy, charge and magnetic flux of 
the magnetized fluid. Using the conservation laws of the circulation of magnetized flow found by 
Bekenstein and Oron for the ideal magnetohydrodynamic (MHD) fluid, and also for the flow with 
zero conducting current, we show that, for nearby equilibria that conserve the quantities mentioned 
above, the relation SQ = is satisfied. We also discuss a formulation for computing numerical 
solutions of magnetized binary compact objects in equilibrium with emphasis on a first integral of 
the ideal MHD-Euler equation. 
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I. INTRODUCTION 

Recent observations of anomalous X-ray pulsars, or 
soft 7-ray repeaters suggest the existence of neutron stars 
associated with magnetic fields strong enough to affect 
their structures in hydrostationary equilibrium (see, e.g. 
[l[). Such objects have not been found in binary neu- 
tron star systems, but hypothetically strongly magne- 
tized neutron stars or black holes may form binary neu- 
tron star or black hole - neutron star systems. In this ar- 
ticle, we model such magnetized binary compact objects 
in close circular orbits, assuming that the spacetime and 
magnetic fields satisfy a helical symmetry and that the 
stars are in equilibrium. 

The helically symmetric spacetime was introduced by 
Blackburn and Detweiler [2j to model binary compact 
objects in close circular orbits in general relativity. In 
such spacetimes, equal amounts of ingoing and outgoing 
radiation are propagating, and hence these spacetimes do 
not admit flat asymptotics, because the steady radiation 
field carries an infinite amount of energy. Nevertheless, 
it is expected that such a spacetime has an approximate 
asymptotic region up to a certain radius, where gravita- 
tional waves are propagating in a curved background, and 
the energy of radiation docs not dominate in the gravi- 
tational mass of the system. Such a solution, however, 
has not yet been calculated successfully in the regime of 
strong gravity. Analogously to Schild's result in electro- 
magnetism for two oppositely charged point particles Q , 
circular orbits of two point particles have been obtained 
in post-Minkowskian spacetimes [|[. More studies for 
the helically symmetric spacetimes have been reported 
by several authors [5l-fl4| . 

In Q (hereafter FUS), thermodynamic laws for heli- 
cally symmetric perfect fluid spacetimes have been de- 



rived. In the first part of this paper, we extend the results 
of FUS to the magnetized perfect-fluid Einstein-Maxwell 
spacetimes with helical symmetry. As in FUS, we use a 
helical Killing vector k a to define a conserved Noether 
current and associated Noether charge Q [l5l - l2l| . With 
an appropriate choice of the current and a constant of 
the electric potential, the charge Q becomes finite and is 
independent of the 2-surface S on which it is evaluated 
as long as the matter and black holes are enclosed in S. 
We obtain the first law, which relates the change SQ to 
the changes in the baryon mass, entropy, circulation and 
electric current of the fluid, and in the area and electric 
charge of the black holes. Its expression corresponds to 
the mass variation formula for stationary axisymmetric 
spacetimes derived by Carter [22|,[23| (Eq. (l37l) below 1 ). 
Concrete calculations for the variation, SQ, associated 
with the classical action for an Einstein-Maxwell theory 
coupled with a perfect fluid carrying an electric current, 
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are summarized in Appendices and [B] to clarify nota- 
tion and conventions. 

When the late stages of binary inspiral are modeled us- 
ing a sequence of equilibrium solutions of helically sym- 
metric perfect fluid spacetimes (without electromagnetic 
fields), the baryon mass, entropy, and circulation of the 
flow, and the area of each black hole are assumed to be 
held constant (sec e.g. [1, 0-1271). Then, the expres- 
sion of the first law becomes SQ = 0, or for asymptot- 
ically flat systems such as the post-Newtonian, or the 



1 The first law Eq. (137 ft is for generic flows that respect the helical 
symmetry. 
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spatially conformally flat systems, SM = £18 J, as a re- 
sult of the conservations of those quantities (FUS). When 
electromagnetic fields and electric currents are present in 
neutron stars, the circulation of magnetized flow is not 
conserved in general. Hence, it is not possible to find a 
sequence of solutions along which the first law is simpli- 
fied as above without further assumptions for the flow. 
In other words, in order to approximate binary inspiral 
just before a merger by a sequence of quasi-equilibrium 
solutions, one needs to introduce a model for the evolu- 
tion of neutron star spins. However, as shown in Sec. lIIIl 
with an electric current introduced by Bekenstein and 
Oron for a class of ideal magnetohydrodynamic (MHD) 
flows ( j28l - [30l |. see also [3l| for non-relativistic magne- 
tized flow), a generalized circulation of magnetized flow 
is found to be conserved. Applying this law - the gener- 
alized Kelvin theorem for ideal MHD - we show that the 
relation SQ = is satisfied along a sequence of helically 
symmetric equilibria of magnetized binary systems, and 
that the relation 5M = 0<5 J holds for asymptotically flat 
systems. 

The above first law can be applied to actual sequences 
of solutions, and hence in the second part of the paper, in 
Sec. II VI andlVl formulations for computing such equilib- 
rium solutions of magnetized binary compact objects are 
discussed. In particular, we discuss the first integral of 
the MHD-Euler equation, which is a key to compute equi- 
libria of neutron stars numerically. Bekenstein and Oron 
[29j have found a first integral of the relativistic MHD- 
Euler equation using the same current for the case with 
ideal MHD irrotational flow, and also for the case with 
the purely convection current. As irrotational flow is con- 
sidered to be more realistic in the final inspiral stage of 
the binary neutron stars and the black hole - neutron star 
binaries [321 ] , we introduce the first integral by Bekenstein 
and Oron for ideal MHD irrotational flow, then derive 
a somewhat different first integral, which may be valid 
only on an initial hypersurface E t , and write down a set 
of equations for the magnetized irrotational flow suitable 
for numerical computations of binary neutron stars and 
black hole-neutron star binaries in equilibrium. 

We follow the conventions and notation in FUS. For 
a one-form w a , the exterior derivative (dw) a p (within 
index notation) is defined by 



(dw) a p := \7 a Wj3 - V/jW a , (2) 



and for a two-form w a f) = W[ a m by 



{dw) al3l := 3V[ a ^ 7 ] = V Q w ( 3 7 +V (3 w 7Q -|-V 7 w Q!/ 3. (3) 



II. THERMODYNAMIC LAWS FOR 
EINSTEIN-MAXWELL SPACETIME WITH 
CHARGED AND MAGNETIZED PERFECT 
FLUID 

A. Zeroth law and constancy of the electric 
potential on the Killing horizon 

We consider a globally hyperbolic spacetime (M,g a p) 
and a vector field k a transverse to each Cauchy surface 
(but not necessarily everywhere timelike). This vector 
generates the one-parameter family of diffcomorphisms 
Xt ■ The action of Xt to a spacclike sphere S on a Cauchy 
surface generates a timclikc surface, T(<5>) = UtXt(S), 
called the history of S. Then, as in FUS, k a is called 
a helical vector if there is a smallest T > for which P 
and Xt{P) are timelike separated for every point P out- 
side of the history T(<S) • Very often, k a can be written 
k a = t a + nct) a , where tt = 2-k/T, t a is a timelike vector 
and (f> a a spacelike vector that has circular orbits with a 
parameter length 2-7T (see, FUS). 

Each Cauchy surface of the helically symmetric space- 
time does not admit flat asymptotics because the energy 
of the radiation generated by a binary equilibrium even- 
tually dominates and causes a divergence. Therefore, as 
discussed in FUS, the future (past) horizon H ± in he- 
lically symmetric spacetime is defined by the boundary 
of the future (past) domain of outer communication , 
where P £ M. is in T) ± if the future (past) timelike curve 
c(A) through P(:= c(0)) remains outside of T(«S) of each 
sphere S for a certain Ao, A > Ao. It is also shown that, if 
the history T(S) of a sphere S is in I? ± , the future (past) 
horizon agrees with the chronological past (future) of the 
history T, H± = dF(T). 

The conditions of the theorems by Friedrich, Racz, and 
Wald [33[ are modified to make them suitable for helically 
symmetric spacetimes. With the assumption that the 
null energy condition holds: R a pl a l^ > for any null 
vector l a , those theorems yield the following properties: 
are Killing horizons, the shear o~ a p and the expansion 
6 of a null congruence vanish on H ± , the Killing vector 
k a is parallel to the null generators of the horizons, and 
the surface gravity k of each disconnected horizon defined 
by 

k p Vpk a = nk a (4) 

is constant on each connected component of H (FUS). 
The Raychaudhuri equation, 

|[ = -R ap rf-2a a ^-\e\ (5) 

is used to demonstrate the above properties. It implies 
R a pl a lP = on the Killing horizons % ± . Assuming there 
exists no material flow through the horizon but there 
exists an electromagnetic field F a p := (dA) a p = V ' a A$ — 
\7pA a , where A a is the electromagnetic potential one- 
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form, we have 

R a pk a kP = 8irT^k a k f3 = 2F a - ( F^k a kP 

= \(E a E a + B a B a ) = (6) 

on 'H ± , where T^a is the stress-energy tensor for the elec- 
tromagnetic field, and E a and B a are the electric and 
magnetic components with respect to the helical vector 
defined by 2 

E a := F aP kP, B a := ^e aM5 F^k s . (7) 

Note that, as a consequence of ©, E a and B a are both 
null on 'H ± . Using the Cartan identity, 



k p {dA) 0a = £ k A a - V a (fc%), 



(8) 



and assuming that A a respects the symmetry £kA a = 0, 
one can introduce an electric potential in the rotating 
frame E a = -V Q $ E . 3 Since E a k a = B a k a = and 
E a and B a are both null on rl ± , E a and B a are nec- 
essarily parallel to the null generator on ri ± . Then, for 
any vector if 1 tangent to 'H ± , r\ a E a = — r] a \7 a $ E = 0, 
which implies that $ E is constant on the Killing horizon 

■h ± mil. 

The potential $ E is defined globally if the domain of 
outer communications is simply connected, and $ E is 
unique up to the constant of integration. The constant 
may be chosen <E> E — > as r — > oo for asymptotically flat 
systems. For the helically symmetric system, we set the 
constant by the condition 



^- j> k^A 1 F af3 dS af) 







(9) 



on the boundary sphere S which encloses all black holes 
and neutron stars, and on which a family of Noether 
charges is defined in the next section 4 . The total electric 
charge of the system is defined by the surface integral over 
the sphere S, 
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-in 



Q K := — <b F af3 dS a p, 



(10) 



2 If k a would be normalized by k a k a = — 1, E a and B a could be 
interpreted physically as the electric and magnetic fields mea- 
sured by the observer of 4- velocity k a . Note however that in 
general k a k a ^ —1; even k a k a = on H^. 

3 One can avoid the assumption that the field A a respects the heli- 
cal symmetry. Eq. Q implies (dE) a p = — fc 7 ((iF) 7a(3 — £fc-F a ^ = 
for (dF) a fj 1 = and the symmetry £fc-F Q) 3 = 0. Hence, from 
the Poincare lemma, 3 <I> E such that E a = — V ce 3 >E if the domain 
is connected and simply connected. 

4 For an asymptotically flat spacetime, the Noether charge defined 
on S with the choice of Eq. (J9J, and then the radius of S taken 
to be r —> oo, agrees with a choice <t> E -> at r — > oo (see, 
Scc. lHBt . 



and the condition ([9]) is rewritten for —k a A a = $ E + C 
with 



C 



(11) 



B. First law for systems with a single Killing vector 

1. Definition of the Noether charge Q 

Given a 1-parameter family of magnetized perfect-fluid 
Einstein-Maxwell spacetimes specified by 

S(A) := lg a p(\),u a (\),p(\),s(\),A a (\) 7 r(\)}, (12) 

a family of Noether charges is defined on any sphere S 
that encloses all black holes and neutron stars associated 
with the electric charge and current [l6r - [2"l| : 



where 
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8^ 



^dS aPl (13) 



Vk 13 + k a< 3 13 (14) 



and 53" (A) is any family of vector fields that satisfies 
1 d 



-gdX 



(15) 



9 Q being defined by Eq. (fA~30]) in Appendix [£] The 
vector 03" (A) is written, 



<8"(A) 



16tt 



(5 q V 5 -5 q V 5 )Ia=o 3 7 *(A) 



+ — F^| A=0 [^(A)-^(0)J +0(A 2 ), (16) 

o 

where is the covariant derivative of the metric g a p{Q) 
and b is a fixed parameter. 

We choose 03" (A) to make Q(A) finite; and, as we will 
see below, Q(X) is independent of the sphere S, as long as 
S encloses the fluid and black holes associated with elec- 
tric charge and current. We first choose the parameter 
b in definition (fTB"|) to have Q(0) satisfy these proper- 
ties. Regardless of the choice of 03" (0), the variation of 
the Noether charge 8Q is finite and independent of the 
sphere 5". The change in the Noether charge SQ results 
in the first law for the Einstein-Maxwell spacetimes with 
charged and magnetized perfect fluid and Killing hori- 
zons, associated with a single Killing vector to impose 
the stationarity of the system. 

In the calculation of the variation SQ, the Eulerian 
change of each quantity in Eq. (fT2"j) is defined by SQ := 
^Q(A), and the Lagrangian change at A = is given by 



AQ=(5 + £ i )Q 1 



(17) 



where £" is a Lagrangian displacement. The definition 
of Lagrangian perturbations is given in Appendix IA 11 
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2. Independence of Q(0) on the location of S 

When the electromagnetic field satisfies F a pF a P = 
in the region where the sphere S is located, b = 1 is cho- 
sen in Eq. (|16|) to make Q(0) finite and independent of 
S. In this case, we have *8 a (0) = 0. When the steady 
electromagnetic radiation is propagating everywhere in 
the spacetime, b = 1/2 is chosen. Then, *8 a (0) be- 
comes <B Q (0) = J F l/3a | A=0 ^ / 3(0)/87T. For the former case, 
a contribution from the gravitational radiation field to 
the charge Q(0) is subtracted, and for the latter case, 
contributions from the gravitational and electromagnetic 
radiation fields to the charge Q(0) are subtracted; Q(0) 
is finite and independent of S as long as it contains the 
fluid and all black holes in both cases. 

To prove that the charge Q = Q(0) is independent of 
the sphere S, we write Q = Qk + Ql, where Qk is the 
Komar charge and Qi an additional contribution related 
to the surface term of the Lagrangian, with 



8?r Js 

Q L = j (k a <&P - k 0f B a )dS, 



(18) 
(19) 



and rewrite Q in terms of integrals over a spacelike hy- 
persurfacc £ transverse to k a . The boundary of E, 



<9£ = SUt Bi, 



(20) 



is the union of the sphere S and black hole boundaries 
Bi, which is the ith connected component of £ n H + . 
Correspondingly, surface integrals over the ith black hole 
horizon Bi are written, 



Qxi = ~ i V a k f3 dS a(3 , (21) 
Q Ll = I (k a <B? -k^ a )dS a0 . (22) 

JBi 

Then, from the identity 

VpV a k p =R a pk p , (23) 

we have 

Qk — ^ Qi<i 

i 

= ~ [ V a k?dS aP = f R a p k?dS a 

= ~h S^ GapkPdSa ~ ~ikJ s RkadSon (24) 

where the integral over the boundary dT, is related to the 
surface integrals with the orientations, J dJ , Q a/3 dS aj 3 = 
(/ s - E §Bi )Q aP dS a p. If F a0 F af} = is satisfied in the 

i 

neighborhood and outside of the sphere S, the vacuum 
Einstein equation is satisfied in the same region. From 



Eq. ((24]) and the choice 53" (0) = 0, Q is then independent 
of the location of S. For the case F a pF a ^ ^ 0, using 

Ql — ^2 Qn 

i 

= [ V /3 (fc Qs 8 /3 ~kP<B a )dS a = [ V^k a dS a 

= L{h v ^ Aa -^ Fa " Fa ") kldS ^ (25) 



have 



i 

~ f (G a p - 87rT F Q ,3)fc^^ - -L f Rk a dS a . 
8tt J s h 16tt J s 



47T 



E i i k^F^ds, 



(26) 



where Tp a p, the stress-energy tensor of the electromag- 
netic field, is defined by Eq. (|A6|) . To derive Eq. (|2l)|) . we 
have used the Cartan identity (|5]l. the symmetry relation 
£ k A a = 0, and Eq. ©. From Eq. ([25]). it is obvious 
that Q does not depend on the sphere S as long as it 
encloses all black holes and neutron stars; all integrands 
of the volume integrals over E in Eq. (|26j) are zero in the 
region where there are no matters and currents, where 
the sphere S is placed. This argument may be clearer by 
using an expression of the Komar charge associated with 
the Lagrangian, Eq. (|B8|) . given in Appendix iBl 



3. First law for the charge Q 

The generalized first law will be obtained by evaluating 
the variation 5Q in the Noether charge in terms of per- 
turbations of the baryon mass, entropy, circulation and 
electric current of each fluid element, and the surface ar- 
eas and charges of the black holes. To find the change 
5Q, we first compute the difference, 



Q-J2® 



(27) 



between the charge on the sphere S and the sum of the 
charges on the black holes Bi. The calculation is per- 
formed in Appendix [B] and results in Eq. (|B15|) . In 
computing the difference (j27)) . we choose two kinds of 
gauge: the first one is to choose Sk a = using the dif- 
feomorphism gauge freedom, and the second one £ = 
for the Lagrangian displacement as a result of the triv- 
ial displacement (see Appendix IBl and FUS). For a per- 
fect fluid spacetime, it has been shown in FUS that the 
quantity (|27[) is invariant under gauge transformations 
that respect the Killing symmetry. For the case with an 
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electromagnetic field, the same invariance under gauge 
transformations associated with diffcomorphisms and the 
U(l) gauge symmetry is shown to hold for the charge Q 
with a contribution from the electromagnetic fields, as is 
discussed below. 

In the black- hole charges Qi = Qki + Qii, Qki is 
calculated in FUS: 



>,Ki 



-!- / V a k fi dS a p = ^-mA 



8- 



B, 



8tt 



(28) 



where Ai is the area of the ith black hole. The Qn is 
made of contributions from the geometry, electric charge, 
and electromagnetic field. The former has been evaluated 
in FUS following dJ: 



SQ Li = I (k a Q p ~ k?O a )dS a p 

JBi 

= -J-SmAi + — I k a F^ a 6A p dA. (29) 
8tt 4tt J m 

For the latter contribution, since k a Fp a = Ep is parallel 
to the null generator kp on we have 

k a F 0a 6ApdA = k a Fp a g^8A^dA 

= k a Fp a {-k p xO -n p k r )SA 7 dA 
= k^F^iWA^dA 



5{k' y A y )F a ^^(k a np - kpn a )dA 



SiWAAF^dSct 



(30) 



where n Q is the unique null vector field orthogonal to Bi 
satisfying n a k a = —1, and 8k 1 = is used. Hence 

8Q Ll = -^-SKiAi + -!- / 8{WA r )F afi dS a p (31) 

87T 47T J Bi 



The contributions from the horizon are Eqs. (|2"5|) and 
3T|) - and the surface integral in the r.h.s. of Eq. (|B15j) . 



J2 ^ 5 j> k~'A r F af> dS a i 



(32) 



Hence the sum of Eqs. (|3T|), ([32]) and the perturbed ([28 
becomes 



47r J Bi 



Sit 



-k,5Ai 



'Bi 
1 



4tt 



—KiSAi + SfSQf, 
8tt 



(33) 



where the total electric charge of the system Eq. (|10j) is 
rewritten using Stokes' theorem: 



Q E = [ j a dS a + V / F^dS af 

JY, ; 47r JBi 



(34) 



and the electric charge on each black hole is defined by 
Of := i- j> F^dScp. (35) 

Note that $f is defined on each Bi by 

= -A a k a = $ E + C (36) 

and is constant. 

Finally, when Einstein's equation, Maxwell's equa- 
tions, their linear perturbations and the equation of mo- 
tion are all satisfied, the first law, which relates the 
change of the Noether charge to changes in the ther- 
modynamic and hydrodynamic equilibrium of matter, in 
the electric current and electromagnetic field, and in the 
area and electric charge of the horizon, is derived from 
Eqs. ([515)1 and 

SQ = [ \- t A(s P u a dS a ) + — ^-A(pu a dS a ) 
is I u u 

+ v l3 A(hup P u a dS a ) - Apk fj A(j a dS a ) 
_ [fk p -fk a )AApdS a 

+ Ef^- 4 ' + ®f S Qf)- ( 37 ) 

Here T is the temperature, s the entropy per baryon, h 
the relativistic enthalpy and v a is defined by the following 
decomposition of the fluid 4-velocity with respect to the 
helical vector: 



u t (k a +v a ) with v a V a t = 0. 



(38) 



Note that u* = u a V a t. In the special case of stationary 
and axisymmetric spacctimes (for which k a is a linear 
combination of the stationary Killing vector and the ax- 
isymmetric one), Eq. ([37p reduces to the mass variation 
formula derived by Carter [13, [23| . 5 

As mentioned earlier, we can verify now that Q{\) is 
independent of the location of the 2-surface S on which 
it is evaluated. In Sec. lIIB^l the charge Q(X) at A = is 
shown to be independent of S, and the variation formula 
Eq. ([37| imply that dQ/dX = SQ is independent of S as 
long as it encloses the fluid and black holes. 



5 To derive the mass variation formula for stationary and axisym- 
metric spacetimes from Eq. I|37|l . one can replace the helical vec- 
tor k a by the timelike killing vector t a . All calculations above 
are valid with this replacement, and now 8Q becomes SM as 
the sphere S goes to infinity (see Sec. Ill B~5l and FUS). Extra 
terms relating to the angular momentum of the fluid and black 
hole (geometry and electromagnetic field) appear in the r.h.s of 
Eq. J3TD as a result of this replacement. 
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4- Gauge invariance of 5Q 

As shown in FUS, for perfect-fluid spacetimes, the 
difference 5(Q — ^2Qi) is invariant under any gauge 

i 

transformation associated with diffcomorphisms that re- 
spects the symmetry k a ] in fact S(Qk — ^2Qki) and 

i 

S(Ql — ^2Qli) are separately invariant Because of 

i 

the contribution from the electric potential at the hori- 
zon, 5(Q — Qi) i s no longer invariant when an elec- 

i 

tromagnetic field is present, neither is each contribution. 
We find, however, that an expression in which the con- 
tribution of the electric charge times the potential at the 
boundary is subtracted, 

8(Q-Y,Qi- ^ J (m k ' yA r Faf>dSa t > )' (39) 

is invariant under the gauge transformation that respects 
the symmetry and the U(l) gauge transformation as 
shown below. 

The gauge transformation associated with a vector 
field r] a is given by 



S n Q = £„<2, f(t,) 



-v 



(40) 



and the corresponding Lagrangian variation is identically 
zero, 



A„ = Sr, + £_„ = 0. 



(41) 

We decompose the vector r\ a with respect to the symme- 
try k a , 

Tj a = r] a V a t k a + f] a , (42) 

with f) a V a t = 0. 

Then, the change in 5(Ql — T^Qhi) becomes 



= / Vge p k a dS a = I S„Cd 3 x 



5 n [Q L -J2Q 

i 

= f \/ a (Cf, a )d 3 x = --L / F 7S F^ S k a f dS a0 , 

(43) 

where we used the relation S V C = £ v jC = V Q (£f) Q ) = 
\7 fj (Ck a V a tfj^), with k a \7 a t = 1. The non-zero contri- 
bution to the Lagrangian density C at the boundary <9£ 
is that of the electromagnetic field C-p- 

Similarly 5(Qk — Qki) is calculated from Eq. (124ft : 



S v (Qk - J2 Q«i) = - ^ S n Ra P kf>dS <* 
°t Joy, 

= --L / [£ k As-W s (k^A 7 )}F a5 fdS a p 
+ T- I F, lS F~< s k a ^ dS a p, (44) 

S7r J SB 



where we have substituted R a p — 87rTp a /3 at <9£ and 
Eq. (|A6p . before using Eq. ©. Finally, the last term 
in Eq. ([39]l becomes 



-U„ [ ^A 7 F al3 dS a0 

47r JdT, 

^S n J V 5 (fc^ 7 F« 5 )dS Q 



2vr 



V s (WA 7 F aS )fdS a p. (45) 



Adding Eqs. (J44J), (j43j), and (gSJl, and imposing £ k A a = 
and V pF a P = at 9S, all terms cancel out : 

S v (Q -Y,Qi- ^ j^WA^F^dSnp) = 0. (46) 

Hence the difference (j3"9")l is invariant under a gauge trans- 
formation that respects the symmetry. 

For the U(l) gauge transformation, we directly show, 
instead of Eq. (|3T)f , the invariance of the difference eval- 
uated at the surface S, 



2 - -!- cf> WA^F^dt 



(47) 



under the transformation with a gauge potential /, 

S f A a = V Q /. (48) 
The change in charge Q with this transformation is 



5 f Q = 6 f Q L = i- j{k a F^ - k p r"")S 1 A 1 dS aP . 



Then, the difference ([4"T]) vanishes 

Sf(Q ~ — I k^A 7 F af) dS, 



(49) 



a/9 



4~ 



k [a F f) ~ t] V 1 fdS af} = 0, (50) 



because integration by part of the r.h.s. of the first equal- 
ity becomes an itegration of a divergence over S that 
vanishes, and an integral of 

3V 7 (/c [Q F' 37 ]) = 2/>V 7J F< 9h + V 7 fc 7 .F Q ' 3 + £ k F al3 

(51) 



0. 



when the Maxwell's equation is satisfied on S where the 
current is zero, and k a the Killing vector. 



5. Asymptotically flat systems 

FUS have derived the first law in a Hamiltonian frame- 
work, and shown the relations between Qk and SQl and 
the asymptotic quantities, the ADM mass M . the Komar 
mass Mk associated with the timclikc asymptotic Killing 
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vector t a , and the angular momentum J associated with 
the asymptotic rotational Killing vector (f> a . In the pres- 
ence of an electromagnetic field, the only difference with 
FUS is the following term in SQl 

I (k a Q§ - k?Q%)dS a p, (52) 

J oo 

where 

(f> := lim f , (53) 

JoO r ^°° JS r 

with S r is a sphere of a radius r, and 0p is the sur- 
face term of the variation of electromagnetic Lagrangian 
defined by 

Q « = J_ F P<*SA . (54) 

However, this does not contribute to SQl, because, for 
asymptotically flat systems, the components of A a are 
0(r _1 ) or lower, and, accordingly, those of F al3 are 
0(r~ 2 ) or lower, hence 

/ (k a @P - k p Q^)dS aP = lim I 9£V Q rr 2 df7 = 

(55) 

where the relations k a V a t = 1 and k a V ' a r = have been 
used. Therefore, as in FUS, 

Qk = -^-l V a k?dS aP = ]-M K -nj (56) 

87T 2 

SQl = I (k a @ p - k p <d a )dS af! 

J OO 

= SM - ^SM K + SQJ (57) 
which results in 

SQ = SAI - nSJ. (58) 

As we will see below, when two nearby equilibria are 
compared conserving the integral quantities, including 
the generalized Kelvin circulation for magnetized flow, 
and the areas and electric charges of the black holes, so 
that the r.h.s. of Eq. (|37|) vanishes, the first law is simply 
written SQ = 0, or SM = O <5 J for asymptotically fiat 
systems. 

III. COMPARING STATIONARY SYSTEMS 
A. Ideal MHD flow 

1. Conservation of rest mass, entropy and electric charge 

For an isentropic fluid, conservation of rest mass and 
entropy can be expressed by means of a Lie derivative 
along the fluid 4-velocity u a : 

£«0\/=<7) = 0, £ u s = (59) 



and if these quantities arc conserved in the perturbed 
states, the perturbed conservation laws have first inte- 
grals 

A(p v / ^) = 0, As = 0, (60) 

Since we assume that the electric current is not necessar- 
ily colinear to the fluid 4- velocity, conservation of electric 
current, 

£,V=5 = 0, (61) 

does not imply another perturbed conservation law anal- 
ogous to Eq. (|60| with respect to the lagrange perturba- 
tion of the congruence of flow lines, that is, A(j a dS a ) ^ 
0. However, its volume integral over the neutron star 
should vanish because of the conservation of electric 
charge: 

SQfn = S f j a dS a = f A(j a dS a ) = 0. (62) 

2. Conservation of magnetic flux for ideal MHD 

Assuming perfect conductivity for the magnetized flow 
of the neutron star matter, the ideal MHD condition 



F a pu 15 = 0, (63) 
is satisfied, and the curl of Eq. (|63| becomes 

£ u F aP = (64) 

as a result of the Cartan identity and (dF) a p 7 = 0. 



Eq. (|64p implies the well known conservation law of mag- 
netic flux, Alfven's theorem. Let us introduce the map ip T 
as the family of diffcomorphisms generated by u a , namely 
the curve r — > i/Jt(P) has the tangent vector u a (P) at a 
point P. For any closed curve c contractable to a point, 
we consider the closed curve c T = ip T oc obtained by mov- 
ing each point of c during the proper time r along the 
fluid trajectory through that point. Then the conserva- 
tion of magnetic flux, which is the integral form of the 
law (|6~4"|) , amounts to the conservation of the integral of 
the 1-form A a along the closed curve c T in the fluid : 

I A a dl a = C F = const. (65) 

J C T 

When the perturbed state also satisfies the ideal MHD 
condition, the perturbed version of the conservation of 
magnetic flux (fT2")) has a first integral 

AF a/3 = (66) 

and hence AF a p = (dAA) a p = is satisfied for any re- 
gion in the fluid. For (dAA) a/ 3 = to be satisfied, it 
suffices that AA a = V a $ for some scalar field '5. Con- 
versely, as long as the fluid support (neutron star) is star- 
convex, the Poincare lemma guarantees the existence of 
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ty. As a result, the last term of the volume integral of 
Eq. ([57)1 vanishes: 

/ -(j a k^-fk a )AA p dS a 

= I -(j a k^-fk a )^dS aP 

+ / Ve(j a kP-jPk a )VdS a = 0, (67) 

because there is no electric current on the boundary sur- 
face d£, and a relation, V /3 (j a fc /3 - j f3 k a ) = £ k j a + 
j a \7 pk 13 — fc Q V pj 13 — 0, is satisfied for the conserved cur- 
rent j a that respects the symmetry. 

3. Conservation of circulation for ideal MHD: generalized 
Kelvin's Theorem 

When two equilibria of some ideal MHD flow are com- 
pared with the same rest mass, same entropy and same 
magnetic flux, the perturbed conservation laws (|60|) and 
(|6"6"]) . as well as Eq. (j6"7| are satisfied. Then the change 
in the Noether charge {3"7} becomes 

SQ = ( [v p A(huppu a dS a ) - Apk A(j a dS a )] 

+ E(^*^ + *w)- ( 68 ) 

For some perfect fluid without magnetic field, the cir- 
culation of the flow is conserved. The curl of the rela- 
tivistic Euler equation u^uip a — is written £, u Cj a p = 
where uj a p is the relativistic vorticity defined by Q a p = 
(d(hu)) Q f} and a corresponding integral law, known as 
Kelvin's theorem, is the conservation of circulation, the 
line integral of hu a along an arbitrary closed curve co- 
moving with the fluid. As shown in FUS, the contribution 
from the circulation to the change in the Noether charge 
5Q is included in the term 

/ v f3 A(huppu a dS a ), (69) 

which vanishes when the circulation is conserved in the 
perturbed flow, for example, when the irrotational flow, 
or the corotational flow, is maintained. This can be 
shown in the same way as eliminating a term (|67p us- 
ing the conservation of magnetic flux. 

The integral in Eq. (|68p . however, does not in general 
vanish for magnetized flows, or even for ideal MHD flows, 
because of the lack of a conservation of circulation law in 
the magnetized case. This can be seen from the relativis- 
tic MHD-Euler equation which is not the inner product 
of the fluid 4- velocity and an exact two form, because of 
the Lorenz force on the right hand side, 

u p (d(hu))p a = -F a pf. (70) 



Nevertheless, Bekenstein and Oron [29[ (see also (30j ) 
have found that, if the 4-current takes the form 

f = Vp{pu a qf 3 - pu?q a ), (71) 

where q a is an auxiliary vector field, one can obtain a gen- 
eralized conserved circulation for magnetized flow. This 
4-current is derived from the variation of a Lagrangian in 
which the ideal MHD condition is added as an interaction 
term pq a F a pu" with the Lagrange multiplier pq a . The 
form ([7TJ) manifestly satisfies the electric charge conser- 
vation: V a j a = 0. Note that, for a given 4-current j a , 
one has the degree of freedom to change q a according to 

q a i y q a + Xu a (72) 

for a scalar A without affecting the value of j a . 

Using V a (pu a ) = [cf. Eq. {59])], the 4-current {711) 
can be rewritten 

j a = £ q (pu a ) + pu a \7pqP. (73) 

Substituting the form into the Lorenz force, we get 

- p F a pf = - p F afi Z q {pu p ) = (dr,) aP u^ (74) 

where 1-form rj a is defined by 

V* ■= Fanqt (75) 

and a relation {Cl[) from Appendix[C] which is implied by 
the ideal MHD condition {63]), is used. Note that, thanks 
to {63]) . the 1-form r\ a does not depend on the specific 
choice of q a within the range allowed by ([72]) . By means 
of ([71]) . the MHD-Eulcr equation ([70]) is simply written, 

u (dw) 0a = 0, (76) 

where w a is the generalized momentum 1-form defined 
by 

w a := hu a + r\ a . (77) 

From Eq. ([76]) one can easily deduce a generalised con- 
servation of circulation law for ideal MHD flows. Indeed, 
defining the vorticity u) a p of a magnetized flow by 

w Q/ 3 = V a Wf3 - Vpw a = (dw) a p, (78) 

the Cartan identity, combined with Eq. (|76[) and the iden- 
tity dui = d 2 w = 0, yields 

£«w Q(9 = 0. (79) 

By means of the Stokes theorem, this conservation law 
can be put in the following integral form [using the same 
notation as in Eq. (|65]) ]: 




(hu a + r] a )di a = C m = const. (80) 
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This law, which has been first derived by Bckenstein and 
Oron [20|, constitutes a generalization to ideal MHD of 
the relativistic Kelvin's circulation theorem (which cor- 
responds to r] a = 0, see e.g. (HJ) 

One can repeat the same argument as for the magnetic 
flux in the previous section. The perturbation of Eq. ([79)) 
for the magnetized vorticity has first integral 

Aw Q(3 = 0, (81) 

which implies Aw a ^ = (dA(hu+rj)) a fj = 0. The Poincare 
lemma guarantees that a function vff exists on the star- 
convex fluid support such that A{hu a + r) a ) = V Q <I'. 

It is also suggested from Eq. ([76)) that an irrotational 
magnetohydrodynamic flow, co a p = 0, is described by a 
velocity potential $ that satisfies 



hu a + 7] a 
and, in this case, $ = A<£>. 



4- First law for the 



MHD with Bekenstein 
current 



(82) 



Oron 



For ideal MHD flow with the Bekenstein- Oron current 
(I7T|) . the first law of the form Eq. ([6"8")l is further simpli- 
fied when comparing two nearby equilibria that conserve 
the circulation of a magnetized flow, (|80|) . Substituting 
Eq. (|7ip to the second term of the integrand of the vol- 
ume integral in Eq. ([6"8")) . we have 

- ApkP A(j a dS a ) 

= (£ k A 7 - k P F M ) A [{pu a qi - pu^q a )dS a ] 

- V 7 {ApkP A [(pu a q~< - pu~<q a )dS a }} 

= v p A( V0 pu a dS a ) - v p q 1 AF MP u a dS a 

fc /3 F (37 A( / 9U 7 g Q dS a ) 

- V 7 {A^ A {(pu a q~< - P utq a )dS a ] } , (83) 

where the relation AV f3 {f af} dS a ) = V (j A(f a ^dS a ), valid 
for any antisymmetric tensor / a/3 , and the Cartan iden- 
tity ((5J) are used, and the symmetry A 7 = is imposed. 

Since the ideal MHD condition ([63)) is satisfied, terms 
including F a pu^ are discarded. Also a term involving 
F a pAv,P is discarded, because Alt' 3 is colinear to u 13 (see 
Eq. (|AT0])l Substituting Eq. (JS3j) to Eq. <J37J) , the inte- 
gral of the last term of Eq. ([53")) becomes a surface integral 
on <9£ which vanishes. Hence, the first law (|3T|) for the 
Bekenstcin-Oron formulation of ideal MHD is written 



—jU^Fpj A(pu a q~ r dS a ) 



SQ = 



— A{spu a dS a ) 

nit 



h - Ts 



A(pu a dS a ) 



■v p A[{hu + r] & )pu a dS a ] 
VVAF^ pu a dS a 

- (j a kP -fk a )AA p dS a 



? [h Ki5Ai 



(84) 



Introducing the following notation 

(1Mb := pu a dS a , dS :— sdM^,, 
dC a := (hu a + r] a )dAlB, 



(85) 



we further rewrite Eq. 



T 

— AdS 



h-Ts 



-AdM B + v a AdC a 



- v^q'AF^ dM B - {j a k - j f3 k a )AAp dS a j 
SA+^QfV (86) 



When the circulation of magnetized flow is conserved, 
there exists a potential ^ such that A(hu a + rj a ) = V tt $. 
Applying an argument analogous to that for the magnetic 
flux in Sec. IIII A 2| a term for the circulation of magne- 
tized flow in the r.h.s. of Eq. ([54")) vanishes: 



v^Aihup +rj P )pu a dS a 



= [ {pu a v p - pu' 3 v a )V^dS a = (87) 



where v a dS a = is used in the first equality, and the 
last equality is proved in the same way as in Eq. (|67[) be- 
cause of a relation, V p{pu a v^ — pu^v a ) = V ( g(/9u ,9 fc Q — 
pu a k p ) = -£ k (pu a ) - pu a V f3 kP + k a V f3 (p U P) = 0. 
Therefore, the rest mass, entropy, circulation of magne- 
tized flow and magnetic flux are all conserved in the per- 
turbation of ideal MHD flow with the Bekenstein- Oron 
current ([7T)) . namely, Eqs. (|60|) . (|66)) . and (|81f) arc satis- 
fied, the change in the Noether charge ([54")) becomes 



SQ 



E 



KiSAi + $f SQ] 



B. MHD flow without conduction current 

It is expected that the inner core of the neutron star 
may be composed of a mixture of superfluid protons and 
high-energy particles. Such flows are well described by 
an ideal fluid without conduction current but only con- 
vection current: 



J 



peu 



(89) 



where e is the electric charge per baryon mass [28f . 
Conservation of rest mass, V Q (pu a ) = 0, and current, 
V a j a = 0, imply that the specific charge e is conserved 
along fluid flow lines, 



0. 



(90) 
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Substituting the current ((89)) into the first law ([37]) . wc 
have 



SQ = I {^T A( - S dS ^ 
+ v p A[(huf3 + eA f3 )pu a dS a 



h-Ts 



u 



A{pu a dS a ) 
A(epu a dS a ) 



(91) 



and also into the MHD-Eulcr equation (jTOjl . 

u p (d(hu + eA)) Pa + Apu^V a e = 0. (92) 

As shown in [28| , the circulation of the magnetized flow 
defined by 



(hu a + eA a )d£ a . 



(93) 



is conserved only when the closed curve c T is taken along 
a curve of constant specific charge e. If we further assume 
that the charge is distributed initially satisfying 



e = e(A a u a ) 



(94) 



(or merely e = constant in the simplest case) , the curl of 
Eq. ([92]) becomes a law of conservation of circulation for 
magnetized flow, 



£ u (d(hu + eA))p a = 0, 



(95) 



and r is constant for any closed curved c T comoving with 
the flow. Then, with the same argument in Sec. IIII A 41 
when nearby equilibrium solutions having the same value 
of circulation T are compared, the perturbed conserva- 
tion law, 



A{d{hu + eA)) 0a =0, 



(96) 



is satisfied. Hence, with Eq. ([96]) . a perturbation of 
Eq. (|90]), 



Ae = 0, 



(97) 



and conservation of rest mass and entropy (|60[) . the first 
law for a flow without conduction current is also written 
simply as Eq. ([88]) . It should be noted that the condition 
e — constant may not be too restrictive for an application 
such as the superfluid proton component in a neutron star 
interior. 



IV. INTEGRABILITY CONDITION FOR THE 
MHD-EULER EQUATION IN IDEAL MHD 

When the stationarity or helical symmetry is imposed 
explicitly on the (MHD-)Euler equation, it is no longer 
an evolution equation. In usual methods {24T - I261 ] . its nu- 
merical solution is calculated using its first integral - a 



sufficient condition for the stationary or helically sym- 
metric (MHD-)Eulcr equation being satisfied. Therefore, 
finding the first integral is a key, and also a restriction, for 
computing equilibrium solutions considered in Sec. IIIII 

As shown in Sec. IIII A 31 when the Bekenstein-Oron 
4-current (|7ip is introduced, the relativistic MHD-Eulcr 
equation for ideal MHD flows takes the form ([76]). If we 
assume that the generalized momentum (|77j) of the mag- 
netized flow respects the helical symmetry, £kW a = 0, 
then a first integral is immediately derived for corota- 
tional and irrotational flows, in a way fully analogous 
with the non- magnetized case [42[ (see also [3f|): the 
Cartan identity k^oop a = £kW a — S7 a (wpk@) reduces 
to k^uipa = — V a (wpk@) and, for an irrotational flow 
(ujfj a = 0), or for a corotational one [u a colinear to k a 
so that ([75]) implies k^ujf; a — 0], we get the first integral 
w a k a = const. 

However, it turns out that the assumption £kW a = 
is too restrictive when applied to a corotating flow. In 
view of ([77]) and ([75"]) . it would yield the first integral 
w a k a = hu a k a + F a pk a qP = const. Now, the colinearity 



of k a and u a , along with the ideal MHD condition (]63[) , 
implies F a pk@ = 0. Hence the first integral would reduce 
to hu a k a = const, i.e. exactly the same as in the perfect 
fluid case, without any Lorcntz force term. 

In Bekenstein and Oron's theory [29|, HO] , the momen- 
tum w a defined by ([77]) and ([75]) contains the Lagrange 
multiplier q a . Because q a is not a physical quantity, it 
does not necessarily obey the helical symmetry. This has 
been noticed by Bekenstein and Oron, but has not been 
taken into account when the first integral was derived. 
In this section, we briefly review the properties of the 
4-current by Bekenstein and Oron, then derive intcgra- 
bility conditions for the case when q a does not respect 
the symmetry. 

A first integral for an axisymmetric and rigidly rotat- 
ing neutron star has been derived by Bonazzola, Gour- 
goulhon, Salgado, and Marck [13 (hereafter BGSM). In 
Appendix [Dl it is shown that the Bekenstein and Oron 
theory can also accommodate the BGSM formulation if a 
term involving £kq a is kept in the MHD-Eulcr equation. 



A. Bekenstein-Oron 4-current 

From ([75]) . the Bekenstein-Oron 4-current can be ex- 
pressed as 



1 



:£ q {pu a y/=g) (98) 



= -p£ u q a + u a V (pqP). (99) 

Let us recall that j a is invariant under a change of q a 
of the form ([72"]) . Without loss of generality, a condi- 
tion such as q a u a = 0, or q a \7 a t = 0, may be imposed, 
although these are not used below. 

The 4-current must obey the helical symmetry, namely 
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its Lie derivative along k a must vanish: 

£kj a = Vp(pu a £ k q 13 - pu?£ k q a ) = 0, (100) 
where £ k q a ^ 0. Using ([95]) and (|M|). we can write 

£feJ Q = -^£ [fe , g] (^ Q V r 5) (101) 
y/— g 

= -p£ u £ k q a + u a Vp(p£ k qP) = 0, (102) 

where the commutator notation [fe,3] Q = £ k q a is used. 
Note the commutation relation £ k £ u — £ u £ k = £[fc «i = 
0, for u a respects the symmetry. In the above expressions 
for £ k j a , it is noticed that we have again the freedom to 
add to £ k q a a vector proportional to u a 7 as £ k q a 
£ k q a + \u a . 

B. Helically symmetric MHD-Euler equation 

We first rewrite the MHD-Euler equation by isolat- 
ing the Lie derivative along the helical vector k a . Using 
the decomposition (|38|) of the 4-vclocity, the MHD-Euler 
equation (|76p divided by u l is written 

{k? +vP)(dw) Pa = -V a (wpk ) + £ k w a + v {dw)p a 
= 0. (103) 

Since r\ a u a = F a pu a q^ = for ideal MHD, we have 

w a u a = (hu a + n a )u a = -h, (104) 

hence 




Substituting this relation into ()103|1 . we obtain 

v «(^ + + £kw a + v /3 (dw) f3a = 0. (106) 

Since both hu a and F a p respect the helical symmetry, 
we have, given the definition (|77p of w a , 

£ k w a = £ k {hu a + Fapq 13 ) = F a p£ k q . (107) 

Hence Eq. (|106|) becomes 

V a (^+w f3 v^ + F a/3 £ fc / + v?{dw)p a = 0. (108) 

Starting from this form of the MHD-Euler equation, let 
us discuss two cases: the corotational flow and the irro- 
tational one. 

a. Corotational flow: The flow is corotational if the 
fluid 4-velocity is parallel to the Killing vector: u a = 
u l k a . This amounts to setting v a = in the decomposi- 
tion ([55]) of the 4-velocity. Accordingly, Eq. (|108[) reduces 
to 

Va(^) + F aP £ k qP = 0. (109) 



Note that, thanks to (|99|) and the ideal MHD condition 
(|63)l . we have 

F aP £. h gP = ~— t F a pf (110) 
pvr 

in the corotating case. 

b. Irrotational flow: In the Bekenstein and Oron 
ideal MHD theory, the magnetized flow is called irrota- 
tional when the vorticity uj a p = (dw) a p defined by (|78p 
vanishes identically. The MHD-Euler equation ([76]) is 
then always satisfied. Via the Poincare lemma, a flow is 
irrotational if, and only if, there exists (locally) a poten- 
tial $ such that w a = V Q <I>. Since w a v a = u Q V Q (f> = 
and v^{dw)p a = 0, Eq. (|108|> reduces to 

V Q (- + £„*) + F aP £ k qP = 0. (Ill) 

Note that, contrary to the corotating case, the contri- 
bution of the Lorentz force is divided into two terms: 
F a p£ k q and the term involving the potential $. 

C. Integrability conditions 

Under the assumption of helical symmetry without any 
restriction on the fluid flow, the integrability condition 
for Eq. (|108[) is that the last two terms in the left hand 
side be the gradient of a function /, 

F a p£ k q P + vP{dw)p a = V a /. (H2) 

It may also be possible that each term is separately inte- 
grable, that is, with two functions / and g, each term is 
a gradient, 

F aP £ k q p = -£ k qP(dA)p a = V a /, (113) 

and 

V P(dw) 0a = V a5 . (114) 

Therefore, the problem of finding a current with which 
the helically reduced MHD-Euler equation has a first in- 
tegral is replaced by the problem of finding the Lagrange 
multiplier q a that satisfies the above integrability condi- 
tions. As mentioned in p9| . however, the vector q a is not 
a freely specifiable quantity, and hence it is not trivial to 
find such a q a , even for corotating or irrotational flow 
where the v^(dw)p a term vanishes and the integrability 
condition reduces to Eq. (|113|) . 

V. FORMULATIONS FOR MAGNETIZED 
BINARY NEUTRON STARS IN EQUILIBRIUM 

A. Bekenstein and Oron's first integral for 
magnetized irrotational flow 

As mentioned earlier, assuming the current is writ- 
ten as in Eq. (|7T|) . and the vector q a respects the sym- 
metry, the MHD-Euler equation is integrablc for irrota- 
tional flow. Since the canonical momentum w a defined 
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in Eq. ([77]) respects the symmetry, £ k w a — 0, and the 
velocity potential for the magnetized irrotational flow is 
defined by Eq. ([82]) . the first integral is written £/.$ = 
constant (which is equivalent to w a k a = constant), or 
more explicitly, from Eq. (jllip . 

— f +£ v $ = E, (H5) 

where £ is a constant. Assuming a one-parameter 
EOS, we have three solvable equations, the normaliza- 
tion condition for the 4- velocity, the first integral, and the 
rest mass conservation equation, for the three variables 
{h, u*, <&}. The equation for $ is derived in Sec. IV Dl 

B. A first integral for initial data of irrotational 
magnetized binaries 

Since part of our motivation for calculating numerical 
solutions of compact binary systems is to prepare quasi- 
cquilibrium solutions that can be used as initial data sets 
for binary inspiral simulations, we assume that the mul- 
tiplier q a can be specified freely on an initial spacelike 
hypersurface £4. Then, when all fields and matter sat- 
isfy helical symmetry, and the vector £ k q a is, at least 
instantaneously, proportional to the helical killing vec- 
tor, the term F a p£ k qP becomes integrable 

£ k q a = £ k q t k a , (116) 

and the coefficient £ k q t is a function of Apk@. Note 
that the assumption (|116| is valid only for irrotational 
flow; for corotational flow F a pu^ = implies F a p£ k q 13 = 
£ k q t F a pk 13 = 0. From the Cartan identity (|SJ) and 
£ k A a = 0, and the assumption (|116[) . the term (|113[) 
becomes 

-£ k q {dA) 0a = £ k g t W a {A p k p ). (117) 
Hence, for irrotational flow, Eq. (|llll) is rewritten 

V Q (^ + £ t ,$) + £. k ^V a {ApkP) = 0, (118) 

and is integrable if there is a function / such that 

£ k q l = f(A p k^), (119) 

so that 

^ + £„$ + J £ k q t d{Apk 13 ) = £, (120) 

where £ is a constant. 

If a data set on an initial hypersurface respects heli- 
cal symmetry permanently, the current should necessar- 
ily be stationary, £ k j a = 0. Substituting Eq. (|116[) into 
Eq. QTU2P , we have 



where we have used the facts that p, or u a respect the 
symmetry, and a relation V a fc Q = 0. When the inte- 
grability condition (|119|) is satisfied, a coefficient of u a 
in Eq. (TT2Tj) vanishes, £ 2 k q l = £ k f{A a k a ) = 0, and 
hence a sufficient condition for stationarity of the cur- 
rent £ k j a = is that the coefficients of k a in Eq. (fT2T|) 
vanish, 

£ u £ k q l = £ u f(A a k a ) = 0. (122) 

This condition is equivalent to the component of the ideal 
MHD condition along k a , 

k a F apU P = -£ u (A a k a ) = 0, (123) 

and is rewritten, on the fluid support of E t , as 

£ v (A a k a )=0, (124) 

that is, A a k a is constant along the spatial velocity v a 
defined by Eq. However, as mentioned above, there 

is no guarantee that solutions calculated from the q a of 
Eq. (TIT9]) satisfies Eq. (Tmf or ((T2l|) . 

As we choose £ k q a to be parallel to k a in (|116p . we 
may further restrict q a so that q a itself is parallel to k a , 

cf=q t k a . (125) 

We substitute (|125p to the current ([75)1 to derive an ex- 
plicit form for the current j a , 

j a = £ qtk ( P u a ) + (nfVptfkP) 

= - P k a £ u q l + pu a £ k q l . (126) 

For example, 

£ k q t = constant (127) 

satisfies the stationarity of the current (|122[) and 

q* = [at + b<p + f q (x A )]k a (128) 

satisfies Eq. (|127|) . where f g is a function of coordinates 
x A A — 1, 2 orthogonal to fc Q , k a V a x A = 0, and a, b are 
parameters that satisfy 

a + bQ, = 1. (129) 

Remember that t parametrizes the foliation and the 
symmetry vector is normalized as k a V a t — 1, and 4> 
parametrizes circular orbits with parameter length 2tt 
and k a V a <j) = Q. 

C. A model with q a — [at + bcf> + f q (x A )]q a 

We next consider a more general form of q a where 
neither q a nor £ k q a is proportional to k a . Separating 
the dependence on the coordinate associated with the 
k a , we assume the form of the vector q a to be 



£U'" = -pk a £ u £ k q l + pu a £\<f - 0, (121) 



q a = [at + b(j> + f q {x A )]q a , (130) 
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where q a respects the symmetry 

£ fe <T = o, 

and hence the relation 



(131) 



(132) 



holds. 

For corotational or irrotational flows, the integrability 
condition (|113|) is rewritten with the requirement that 
there exists a function / such that 



F a pf = -q P (dA)/3 c 

or using the Cartan identity, 



V Q /, 



£qA a = V a {A p q t3 - /). 



(133) 



(134) 



When stationarity is imposed to the current, substituting 
Eq. to Eq. ((100)1 . we have 

£kj a = Vp(pu a f - = 0. (135) 

Then, from Eq. (|130|) and Eq. ((7T1) . the current becomes 

r = - puV)V / ,[a* + ty + / g (:c j1 )]. (136) 

c. Corotating flow: This model can be applied to 
corotating flow, as long as one can find a particular form 
of q a that satisfies Eq. (|133[) as well as the stationarity 
and ideal MHD conditions consistently. For corotating 
flow, u a — i/k", Eq. (|136[) becomes 



j a = pu t k a q^p[at + bd> + f q (x A ) 



pu l q a . 



(137) 



Assuming f q (x A ) = and using a + D,b = 1, the combi- 
nation of t and 4> components — Oj* becomes 

j*-fij* = pu^^f - k P q a )W a (j)Wpt 

= -pu 1 ^ - net). (138) 

As discussed in Appendix [Dj when the system is station- 
ary and axisymmetric, and if q a satisfies 



(139) 



the formulation becomes the same as that of [33| for a 
magnetized rotating neutron star. 

d. A trivial model for the irrotational flow: When 
q a is taken to be parallel to k a , 



q f k a , 



(140) 



with £fc<? = 0, the first integral is derived as in the previ- 
ous section, if q l is a function of A a k a ; Eq. (|1 33|) becomes 

- f(dA)p a = ^V a (A p k p ) = VJ. (141) 

The current (|136|) in this case is written 

j a = (0M t K-fc Q ^V /3 [^ + / (? (.T A )]}g*. (142) 

A trivial solution to the condition (|14ip is 



constant. 



(143) 



D. Equation for the velocity potential $ 

To write down an equation for the velocity potential 
$ for magnetized irrotational flow used in an actual nu- 
merical code, we introduce a 3+1 decomposition of the 
spacetime. In this section, spatial indices are Latin. The 
spacetime A4 = K x X is foliated by a family of spacelike 
hypcrsurfaccs (£ t ) te R parametrized by t. The future- 
pointing unit normal to the hypersurface S t is defined 
by n a = — aV Q i, where a is the lapse function. Then 
the generator of time translations in an inertial frame t a , 
and rotating frame (helical vector) k a are related to n a 
by t a = an a + /3 a and k a = an a +ui a respectively, where 
/3 a and u a denote a spatial shift vector in each frame, and 
are related by uj a = (3 a + Q<f) a . The spatial metric "f a b{t) 
induced on E t by the spacetime metric g a p is equal to the 
projection tensor orthogonal to n a , jap = 9a/3 + n a np, 
restricted to St. In a chart (t,x a ), the metric g a p has 
the form 



ds" 



a 2 dt 2 + lab {dx a + I3 a dt)(dx b + P b dt). (144) 



The covariant derivative associated with the spatial met- 
ric 7ah is denoted by D a . 

In the formulation for irrotational flow, the number of 
independent variables becomes three [2{| |42j ■ As inde- 
pendent variables, we choose the rclativistic enthalpy per 
baryon mass, the time component of the 4-velocity, and 
the velocity potential, {h, u*,<f>}. For the first two vari- 
ables, the first integral Eq. (|1 18|) and the normalization 
of the 4-velocity u a u a = — 1 are solved. Using a relation 
derived from Eqs. l(3"8)) and 



(145) 



these equations are rewritten, 



v a D a § + / Wd^fc") = £, (146) 



h 2 [(cm 4 ) 2 - 1] = (D a § - ?f)(D a $ — Tj a ), (147) 

where rj a is a spatial projection of r] ai r\ a = r y a a Va- 

An equation to calculate the velocity potential <3? is 
derived from the rest mass conservation law, Eq. (l59l) . 



— —^(pitay/j) 
-Daiap^v 11 ) = 0. (148) 



Substituting Eq. ()145[) in the above relation, we have an 
elliptic equation for <£, 

D a D a $ = D a ( V a +hu t u a )-(D a <I>-r ]a -hu t uj a ) — D a ^. 

ap n 

(149) 

This equation is solved with a Neumann boundary condi- 
tion to impose the fluid 4- velocity u a to follow the surface 
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of the star. The surface is defined by the vanishing pres- 
sure p = 0, where the relativistic enthalpy is chosen to 
be h = 1 which is always possible when a one-parameter 
equation of state is assumed. Hence, the boundary con- 
dition is written 

u a \7 a h = at h=l. (150) 
and, using E^h = and Eq. (|145p . it is rewritten, 

(O a $ - if - hu*uj a )D a h = 0. (151) 
where V a /i and D a h are normal to the stellar surface. 

VI. DISCUSSION 

A. First law associated with the Bekenstein and 
Oron Lagrangian 

The Lagrangian density of the Bekenstein and Oron 
ideal MHD theory [2^| is based on Schutz's Lagrangian 
density for relativistic fluids [43| . Our Lagrangian density 
for a relativistic fluid C m = —e^/—g (|Al8j) . and the La- 
grangian variation applied to it, is equivalent for the pur- 
pose of deriving the first law. Then, we rewrite the La- 
grangian corresponding to that of Bekenstein and Oron 



C=^—R-e- —F^ + F a0P u^) V~9, 

(152) 

in which the interaction term is replaced by a term F a pu a 
times the Lagrange multiplier pq a which enforces the 
ideal MHD condition F a/3 u a = 0. 

Associating this Lagrangian with the charge Q (fT"3l) . 
we can derive the first law; a calculation of the varia- 
tion 5Q is shown in Appendix [E] Now, the derived first 
law is for the ideal MHD flow, while our first law (|37l) is 
valid for more general MHD flows. Obviously, the argu- 
ment in Sec. Mil applies to the case with the Lagrangian 
(|152[) . Hence, if a sequence of magnetized binary solu- 
tions in equilibrium is constructed assuming conservation 
of rest mass, entropy, magnetized circulation, magnetic 
flux, black hole surface area and charge for a black hole - 
neutron star binary, the first law in the form 5Q = 0, or 
5M = £15 J for asymptotically flat systems, is satisfied as 
for non-magnetized ones, and for the latter case, one can 
apply a turning point theorem to locate a point where 
the stability of solution changes [44 1 . 



B. First integral of MHD-Euler equation 

As mentioned in Sec. lIVl a first integral of the relativis- 
tic MHD-Euler equation is almost crucial for developing 
a successful method to compute equilibrium binary so- 
lutions numerically. When we derive a first integral, we 
need to specify a form of the vector q a , which should be 



consistent with the stationarity as well as the ideal MHD 
condition. However, since q a is not a freely specifiable 
vector, it is not guaranteed that a set of equations ad- 
mit such a q a as solution in general. Also a difficulty to 
have a helically symmetric irrotational binary solution in 
ideal MHD may be explained physically as follows. Be- 
cause the magnetic flux is frozen into the fluid for ideal 
MHD, when the binary system is seen in the rotating 
frame, a poloidal component of the magnetic field may 
be winded up, since the neutron star is spinning in this 
frame. This argument does not rule out the possibility 
to have a helically symmetric magnetized binary neutron 
stars, although it is not trivial at all to find a q a that 
gives such solutions. 

In Sec.|Vl we discuss a formulation for computing equi- 
librium solutions of magnetized binary neutron stars and 
a possible candidate for a first integral of the relativistic 
MHD-Euler equation in ideal MHD flows. Our proposal 
is to assume £kq a be proportional to the helical vector 
k a . It could be possible that this condition is violated 
as the solution is evolved in time, that is, a solution cal- 
culated from the first integral in the Appendix might 
not respect the helical symmetry or the ideal MHD con- 
dition. It would be applicable, however, for computing 
initial data for merger simulations of magnetized com- 
pact objects, because it may be allowed to freely specify 
£kQ a , at least instantaneously on an initial hypersurface. 

In Sec. [V] we also write down a set of equations to be 
solved for an equilibrium of irrotational neutron star in a 
binary system. The formulation for solving the Einstein 
and Maxwell equations are not presented in this paper. 
In usual ideal MHD simulations, the electric current j a 
does not contain dynamical degrees of freedom and, ac- 
cordingly, the Maxwell equation becomes an evolution 
equation for the magnetic flux density. This equation is 
again hard to integrate when the stationarity condition is 
imposed. Therefore our plan is to choose the electromag- 
netic potential one form A a as a variable and to write the 
Maxwell equation as a set of elliptic equations. These el- 
liptic equations can be solved with the same numerical 
method we have developed to solve for the metric poten- 
tials of gravitational fields [HI, HH ■ Our next project is 
to develop such a numerical code. 
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Appendix A: Variation of the Lagrangian 

We begin with a classical action for an Einstein- 
Maxwell theory coupled with a perfect fluid carrying elec- 
tric current, 



s 



Cd x, 



C — £g + C m + Cp + C\ 

16-7T 16-7T 



(Al) 



-g- 

(A2) 



We first define the Lagrange perturbation for the fluid. 

1. Lagrange displacement 

We describe a perfect fluid by its four- velocity u a and 
stress-energy tensor 

T aP = eu a u fs +pq a ^ (A3) 

where p is the fluid's pressure, e its energy density, and 

(A4) 



is the projection tensor orthogonal to u a . We assume 
that the fluid satisfies an equation of state of the form 



P = P{P,s), e = e(p,s), 



(A5) 



with p the baryon-mass density and s the entropy per 
unit baryon mass. (That is, p := mju, with n the num- 
ber density of baryons and nig the mean baryon mass.) 
The electromagnetic stress-energy tensor is given by 



rpa/3 _ _J_ ( pay p/3 _ }_ a(3 p p~/8 



(A6) 



where electromagnetic field 2-form F a p relates to the po- 
tential 1-form by 

F a[j = (gL4) q/3 = \J a A fi - \J fi A a . (A7) 

Given a family of magnetized perfect-fluid Einstein- 
Maxwell spacctimcs specified by 

Q(A) := [g a0 (X),u a (X),p(\), s (X),A a (X) 7 r(X)}, (A8) 

one defines the Eulerian change in each quantity by 5Q := 



We introduce a Lagrangian displacement £ Q in the fol- 
lowing way: Let Q := Q(A), and let be a diffeomor- 
phism mapping each trajectory (worldline) of the initial 
fluid to a corresponding trajectory of the configuration 
Q(A). Then the tangent ^ Q (P) to the path A ^\(P) 
can be regarded as a vector joining the fluid clement at P 
in the configuration Q(X) to a fluid element in a nearby 
configuration. The Lagrangian change in a quantity at 
A = 0, is then given by 



dX 



*-aQ(A)| 



A=0 



(5 + £ e )Q. 



(A9) 



The fact that &\ maps fluid trajectories to fluid tra- 
jectories and the normalization u a u a = —1 imply 



Au a = -u a u^u<Ag M . 



(A10) 



2. Variation of Lagrangian 

Although the variation of the Lagrangian density (|A2[) 
is well known, those calculations are summarized below 
to clarify notation and conventions. A surface term O" 
is kept for the calculation of the first law in Sec. IIIBI 

The variation of the Einstein- Hilbert Lagrangian is 
written as 



1 



-SCr 



1 

167T 



G af3 5g aP + V Q 6g 



_L_ {g <x~ /g PS_ g «P g y6 )V(jdg ^ 



(All) 



(A12) 



The variation of the perfect fluid Lagrangian is de- 
scribed by the Lagrange perturbations. Considering gen- 
eral perturbations in which the entropy and baryon mass 
of each fluid element are not conserved along the family 
(2(A), one obtains 

— = ^u^puo^—g) - \q af) Ag ap] (A13) 

and the local first law of thermodynamics for the fluid, 

Ae = pTAs + hAp, (A14) 
with the relativistic enthalpy h defined by 

P 



h = 



yields 



Ae 



pT 



Hence, we have 



Ae = pTAs 



1 



Ap_ 
P ' 

1 



(A15) 



(A16) 



=hu a A(pu a ^) - -(e+p)q ap Ag aP . 
9 2 

(A17) 
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From these relations, the variation of the Lagrangian 
density for a perfect fluid 



(A18) 



becomes 
1 



-9 



5C V _ 



-Ac - e^g a ^Ag a0 + V a (ef) 



= -pTAs 



-fJ 



-.hu a A(pu a y/—g) 



with the surface term 

0° = (e+j^'V 



(A19) 



(A20) 



TTie variation of the Lagrangian for the electromag- 
netic field 



(A21) 



is calculated as 



1 



1 



l&TTyf—g 
1 



16tt 



2(dSA) aP F^ + 2F ai FfP8g^ 



2F lS F^ l -g^5g a 



- \-Tf5g aP - ^ l3 F a ?5A a + V Q 9£, 



where 0p is defined by 



1 



0° = —F Pa 6A 
An 



(A22) 



(A23) 



TTie variation of the interaction term between matter 
and the electromagnetic field, 



A = A a j a V=9, 



(A24) 



becomes 



jLtfjCi = SA a j a + A a -^A{f^g-) 
<-9 V^g 
l 



Using the relation 
1 



A a -=£s(j a ^). 



wv^g) = v^r^-zn+r^/, (A26) 



(A25) 



-.'/ 



we have 
1 



1 



SCi = j a SA a + A a ^A(f^) 



-g v-9 

+ e [F a pf - A a X7 p f] + W a Q?, (A27) 



where the surface term is defined by 

e? = Ap{re-fn- 



(A28) 



Variation of the Lagrangian density: Finally, the 
above terms are collected and the variation of the La- 
grangian density (|A2|) is derived, 



1 



-.5£ 



1 



■■(SCc 



6£» 



5Ci 



SCi 



l~9 
-pTAs 



-.hu a A(pu a y^g) 



A r 



1 



-9 



Sg a p 



— \G afi -8tt(T q ^ + T" /3 ) 
16tt L v f ' . 

±-(y F a P-4Trr)8A a 

C [V Tj-F aP f + A a Vpf] 
V Q 9 Q , (A29) 



where the surface term <d a is defined by 

& a = eg + + e£ + ef 



g^g^)V P 5g lS + ^ a 5A 



16-7T 



Appendix B: Calculation of S(Q — ~}2Qi 



(A30) 



In calculating a contribution from the volume integral 
to the charge (|27|) . we restrict the gauge in two ways: 
We use the diffeomorphism gauge freedom to set 8k a = 
0. The description of fluid perturbations in terms of a 
Lagrangian displacement £ a has a second kind of gauge 
freedom: a class of trivial displacements, including all 
displacements of the form fu a , yield no Eulerian change 
in the fluid variables. We use this freedom to set At = 0. 
Because 5t = (t is not dynamical), this is equivalent to 
the condition £* = 0. The relation (|A10|) now implies 



1 



=f = =-u a u p Ag a 



(Bl) 



Then, from Eqs. ([AlOf and ([Bljl. we have Au a = 
Au\k a + v a ), while, by u a = u\k a + v a ), Au a = 
A[w*(fc Q +v a )}; thus 



A(k a +v a ) = 0. 



(B2) 
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Then, in the variation of the Lagrangian density (|X29t . 
a term involving a perturbation of the rest mass density 
is rewritten 

hu a A(pu a y/=g) = - Aa0ju*V=5). (B3) 

IT 

To find the change SQ in the Nocther charge, we first 
compute the difference, 



(B4) 



between the charge on the sphere S and the sum of the 
charges on the black holes Bi. 

The difference in the Komar charge Eq. is associ- 
ated with the Lagrangian density as 



Qk — ^2 Qki 

i 

1 



16tt 



R 



1 

167T 



F a pF afs 



f (T a +T F a p)k?dS a 



F a pF af} - A a j a J fc^dS*. 



16tt 

- F- I [G a p- MT a p + T F a p) } k?dS a . (B5) 

Using the relations 

-T a pk p dS a = ~T a B (k f3 + v p )dS a +T a g v dS a 

= ek a dS a + (e+p)u a upv^dS ai (B6) 

and 

-T F a p k^dS a - (^-F aP F^ + A a j°^J VdS 1 
= - ^-F a ~< [£ k A 1 - V 7 (fc^A /3 )] dS a + A a j a k^dS 1 
1 F a ^£ k A 1 dS a + ^-\/ 1 {F a ^k !3 A p )dS a 



in 
1 



4tt 



4^" 7 ' 
Eq. (|B5[) is rewritten 



(B7) 



'if 



fif / 



= - £d 3 x + / {e+p)u a u p v p dS a 



4?r J e 4tt j ss 

+ / Ap{fk a - j a k^)dS a 



WA-F^dS, 



a/3 



^JjG a p - 8 7 r(T Q (3 + r F ^) ] k ff dS a 

±- / fc^ 7 (V,3F^-4^r)^ Q . (B8) 



The variation of Eq. (|B8|) is then 
«(0*r-X)Qxi) 

i 

= - [ 6£d 3 x + I A [{e+p)u a upv f3 dSa\ 

-J-S I F a ^£ k A y dS a + -1-5 [ k^A^F a0 dS a 
4tt 7 s 4tt J as 

+ / A[A (fk a -j a k^)dS a ] 
_ _L S [ [G a p - 8w(T a + Tp a p) } k?dS a , 
-5 [ k~«A 7 (V 8 F al3 ~4nj a )dS a . 



4tt j s 

The integrand of the second term becomes 

A [(c + p)u Q ii0i^dS Q ] 
= A(pu a dS a ) + v fj A(hup)pu a dS, 

+ {e + p)u a u f) £ k f dS a , 



(B9) 



(BIO) 



where Av^ = —Ak 13 = £k^ was used, and the integrand 
of the fifth term is 

A[A fl (jf>k a -j a k f> )dS a ] 
= AAp(fk a -j a kP)dS a 

+ Ap-^A(fV=9)k a dS a 
-Apk A(j a dS a ) 

+ A (j a £ k ^-f£ k e)dS a (Bll) 

where Ak a = £^k a = — £ k £ a , because of our gauge 
choice 5k a = 0. 

The variation of Ql — Qli is given by 



i 

= <f (fc^e* 9 - k e a )dS aP 

JdT, 

= f v ( g(fc Q e /3 -k p e a )ds a 



V 8 @' s k a dS a 



)jo r K'ao a — I £ k Q a dS a , 



(B12) 



where we have used the relation V ' a k a = to obtain the 
last equality. The integrand of the last term in Eq. (IB12|) 
is written as 

£ k O a dS a = {e+p)q a p£ k f dS a 

+ A p {j a £ k e -f£ k i a )dS a 
= (e+p)u a u p £ k ^dS a 

+ A (j a £ k ^ - / £ k C ) dS a , (B13) 

where we used the fact that £ a as well as its Lie deriva- 
tive along k a is spatial £ k £ a V a t = 0. These two terms 
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in Eq. (|B13[) cancel out with the last terms of Eqs. (|B10|) 
and (|B11[) . Note that the current j a respects the sym- 



metry £kj a = 0. 

Finally, we obtain an expression for 8{Q — ^ Qi 



T 



:Aspu a dS a 



A(pu a dS a ) + v p A{hup)pu a dS a 
- A (j k A{j a dS a ) - (j a hP - fk a )AA dS a \ - —8 f F a ^ £ k Ar,dS a + — 5 f WA-F^dSat 

J 47T J S 47T J as 

^-8 [ [G a f! -8ir(T a p +T F a p )}kl 3 dS a - -L 8 [ W A 7 (V - 4*j a )dS a 

OTT J E 47T J E 



1 

i 16tt 



G a/3 - 87r(T a ' 3 + T^) 



4tt 



(V^ - 4nj a )8A a + C [V fj Tj - F a0 f] \ W dS 7 . (B14) 



Note that k a dS a = \/—gd 3 x. When the field equations, their perturbations, and equations of motion are satisfied, 
using £kA a = 0, and Eq. © noting J dY , = <f s — J2i §nv ^l* dEHJ is rewritten 



+ hupyP 



A(pu a dS a ) + v /3 A(hu l3 )pu a dS a 



Apk 13 A(j a dS a ) - {j a kP - fk a )AA dS a \ - V — 8 I k^A 7 F a ^dS a p. (B15) 

J ^ 4?r 7 Bi 



Appendix C: Calculation of Eq. (|74[) 

A relation used in Eq. ([74| is proved in [30j], which 
is repeated here for a reference. Consider a closed 2- 
form ((gLF) Q( 3 7 = 0), and a vector N a such that 
FapN 13 = 0. Then, for any vector q a , a relation 

(dr/) a/3 A' 3 = F a0 £ q N? (CI) 

is satisfied, where rj a is defined by r\ a = i 7 ^ <j^. This can 
be shown as follows: 

(dr)) af) NP = {d{F-q)) aP N p 

= [(q ■ dF) a p - £ q F aP ] 

= F a0 £ q N?. (C2) 

The Cartan identity was used in the second equality and 
the relation F a pN^ = in third one. 



Appendix D: First integral of MHD-Euler equation 
in BGSM formulation 

A formulation for uniformly rotating axisymmetric 
stars with poloidal magnetic fields is derived in [33] ■ In 
this section, we show that the Bekenstein and Oron for- 
mulation of ideal MHD includes a first integral of the 
MHD-Euler equation derived in the BGSM formulation, 
assuming the same symmetry and suitably choosing an 
auxiliary vector q a in the current (|71[) . 



In the BGSM formulation, a stationary, axisymmetric 
and circular spacetime is assumed. And more specifically 
the flow field of rotating star is assumed to be uniform; 
with a constant angular velocity f2, 4- velocity is written 
u a = u l k a = it 4 (f a + VL(f> a ) where t a and 4> a are killing 
vectors. 

Carter has shown [22[ that in stationary, axisymmet- 
ric and circular spacetime, the vector potential and the 
current are such that A a = A t S7 a t + A^7 a (j) and j a = 
j l t a + j^(j> a respectively. Since the vector potential A a 
is assumed to respect the symmetry £kA a = 0, the ideal 
MHD condition F a pu^ = implies, for a corotating flow, 

F a pk? = -£ k A a + V a {A p kf) = V a {A p kP) = 0, (Dl) 
hence 

A a k a = A t + VlA^ = constant. (D2) 
Using this relation, the vector potential is written 

A a = A 4 {V a <f)-nV a t). (D3) 
Note that V a (j>— f2V Q i is orthogonal to the helical vector, 

k a (v a (t)-ny a t) = o. 

Rewriting the current as 

j a =j t k a + Jc/> a , (D4) 

with J = j a (V a <t> - nV a t) = j* - Qj*, the Lorenz force 
becomes 

-F aP f = -F a p<^ = - [-£^A a +\/ a (A^)] . 
p p P L 

(D5) 
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Then, with the symmetry £ ( j > A a = 0, the MHD-Euler 
equation (|T0|) is written 



J 



k (d(hu)) 0ct = —V a {A p ^) 



pu L 



(D6) 



or using 



V Q (-r J - ^tV q A = 0. 
\u l J pu l 

Hence, an integr ability condition of this equation is 
J 



pu L 



(D7) 



(D8) 



Equation (|D7[) and MHD-Eulcr equation for the co- 
moving flow (|109p with the current (|7Tj) agree if the re- 
lation 



F a p£ k q 



J 

pu l 



(D9) 



is satisfied. For example, if the vector q a satisfies 



J 

pvt 



r = fiA^r, (dio) 



the Bekenstein and Oron formulation becomes the BGSM 
formulation [cf. Eq. (|D5j)]. 

Appendix E: Calculation of S(Q — Y^Qi) f° r the 

i 

Lagrangian with Bekenstein and Oron's interaction 
term 

In the Bekenstein and Oron theory, the ideal MHD 
condition F a pu^ = is imposed by adding a constraint 
to the Lagrangian density with a Lagrange multiplier q a , 



A = F a ppu a q 



-<)■ 



(El) 



This term replaces an interaction term, A a j a s/—g, of the 
field and electric current. The variation of C\ becomes, 



1 



-.SCi = -V p 6A a {pu a q f3 - p U Pq a ) 



+ -/=*W> [&{pu a qP^—g) - £ c (pu«VV=fl)] ,(E2) 
The last term is calculated as 

- -^F aP ^{pu a q' 3 ^—g) 

= pu a qP[C(dF) 7a p+d(t;-F) a p} 



V Q (F 01 puiq^ a ), 



(E3) 



where the Cartan identity for the 2-form F a p : £^F a p = 
(,' y (dF) ia i3 + (g?(£ • F)) a p is used, and 

pu a qP(d(t-F)U = ( P u a qP -pueq a )V a (CF^) 
= i a F^ f + V Q [ (pu a qP - p U Pq a )CF lP ] . (E4) 



Hence, using j a = V p(pu a q" — pw q a ), we have 



1 



-.SCi 



1 



:F a ^{pu a qP^—g) 



+ j a 5A a + C [Fcpjf + pu f3 q^(dF) a ^] + W a Qf, 

(E5) 

where 

Of = (pu a qP - pu l3 q a )5A f3 

- (puVf + pu^q-'C + pu'q a ^)F M . (E6) 

To calculate the difference of Noether charge 5(Q — 
YlQi)i we nrs t associate Qk — ^2Qki with the La- 



grangian (|152j) as 

Qk — ^ Qki 



I Cd 3 x 



s 
1 

4-7T 



(e + p)u a u p v p dS a 



F a ~<£ k A y dS 



1 

4tt 



+ / F al 3 P u a q p k~ t dS 1 
is 



WA^F^dS, 

as 



1 [ [G a p - 8ir(T a p +T F a p )}k?dS a 
- / WA^iVpF* -A-Kj a )dS a , 



(E7) 



which corresponds to Eq. (|B8[) . The integrand of the fifth 
term in the r.h.s. of Eq. (|E7[) is rewritten 



F a0 pu a q (3 k' y dS 7 = F a/3 p(k a +v a )q f) v?dS 1 
= k a F olP q p pu'dS 1 + v a g a pu' 1 dS 7 , (E8) 

and combined with the sixth term as 

v a g aP u l3 dSp - W a {k p A )q a pu l dS 1 - k p A p j a dS a 
= v a g a pu p dSf) - £ q {k A p pu'dS 1 ) (E9) 



where Eq. (|98j) and £ k A a = were used. The integral 
of the last term of Eq. (|E9[) over S is rewritten a surface 
integral over 9E that vanishes, because of the gauge in- 
variance under the transformation q a — > q a + Xu a which 
can always be used to make q a spatial, q a V a t = 0. 

The third line of Eq. (|E7j) is replaced by v a g a pu^dSp, 
then a variation of the charge is calculated. 

A difference from the calculation of 5(Qk — Qki) in 

i 

Appendix iBl is the terms, 

A [(e+p)u a u p v p dS a ] + A(v f3 gppu a dS a ) 
= {hup + r]p)v fi A(pu a dS a ) + v A(hup + gf;)pu a dS a 
+ (e+p)u a u p £ k ^dS a + £ k ^r) pu a dS a , (E10) 
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where Av@ = — Afc' 3 = £ k ^ is used. In the calculation 
of 5(Ql — J2Qn), a term £^8° dS a becomes, 

i 

£ k O a dS a = (e+p)u a up£ k ^ dS a 

+ £ k CF^{fm a qe - P uPq a )dS a 
+ (5A P + ^F lP )(pu a £ k q ri - pu {3 £ k q a ) dS a 

(Ell) 

where £ a and £ k £, a are both spatial. The first term and 
a part of the second term in the r.h.s. of Eq. (|E11I) cancel 
out with the last two terms in Eq. (|E10[) . With the Car- 
tan identity, C 1 F 1 p = £^Ap — V^(^ 7 A 7 ) the last term of 



Eq. (jElip becomes 

(6A + CF^)(pu a £ k q f:) - pu^£ k q a ) dS a 
= [AA P - V^(fA 7 ) ] ( P u a £ k qP - pu^£ k q a ) dS a 
= AAp (pu a £ k q f3 - pu?£ k q a ) dS a + C^£ k j a dS a 
- Wp [CA 7 (pu a £ k qP - P uP£ k q a )} dS a , (E12) 

where the second term of the r.h.s. of the last equality 
vanishes for the symmetry, £ k j a = 0, and an integral of 
the last term over S vanishes for the Stokes theorem. 
Finally, we obtain an expression for 6(Q — Qi) f° r 

i 

the Bckcnstcin and Oron ideal MHD theory: 



s(Q-Y,Q*) 



T 

—As pu a dS a 



— + {hup + r\p)v^ 

1 1 " 



A{pu a dS a ) + v p A(hup + ■n P )pu a dS a 



- {pu a £ k q p - pu^£ k q a )AAp dS a \ + 1-8 I If AyF a P dS a p + I F a p 
1-6 f [G a p-8ir(T a p+T F a p)]k l3 dS a -1-6 f WA^pF"? - ^j a )dS a 

OTT Jz 47T J s 



— A{q a pu 1 dS 1 ) + £ k Cpq' < dS 1 



s 
1 

i 16tt 



G ap -8tt(T^ +Tf) 5g a p + -^(V pF a ? - 4nj a )6A c 

47T 



r [ vpTj - F a pf - /vy ] | k 5 ds s . 



(E13) 



This expression is compared with Eq. (|B14I) . Note that, nctized flow explicitly appears as in Eq. 
in the second line of Eq. (|E13|) , the circulation of mag- 
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